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Abstract 

We analyze the radiative leptonic decay mode: B^ — * £2/7 (£ = e, yu) using 
a QCD-inspired constituent quark model. The prediction: B{B(. —>■ iw-f) ~ 
3 X 10~^ confirms that this channel is experimentally promising in view of the 
large number of Be mesons which are expected to be produced at the future 
hadron facilities. 
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Recently, the CDF Collaboration at the Fermilab Tevatron has reported the observa- 
tion of the Be meson, the lowest mass be (be) bound state, through the semileptonic decay 
mode —>■ J/ipi^X The measured mass and lifetime of the meson are 

Mb, = 6.40 ± 0.39 (stat) ±0.13 (syst) GeV (1) 
t{Bc) = OAQtoil (stat) ± 0.03 (syst) ps . (2) 

The particular interest of this observation is related to the fact that the meson ground- 
state with open beauty and charm can decay only weakly, thus providing the rather unique 
opportunity of investigating weak decays in a heavy quarkonium-like system. Moreover, 
studying this meson, important information can be obtained, not only concerning fun- 
damental parameters, such as, for example, the element Vcb of the Cabibbo-Kobayashi- 
Maskawa mixing matrix, but also the strong dynamics responsible of the binding of the 
quarks inside the hadron. Understanding such dynamics is one of the most important 
issues in the analysis of heavy hadrons p|. The observation of the B^ meson at the Teva- 
tron confirms that B^ physics will gain an important role at the future hadron facilities, 
where a large production rate of such particles is expected; in particular, at the Large 
Hadron Collider (LHC), which will be operating at CERN, it is estimated that 4.5 x 10^° 
B^ mesons will be produced per year for a machine luminosity of C = 10^^ cm^^ sec^^ 
at 14 TeV g]. 

Be meson decays can be classified according to the mechanism inducing the processes 
at quark level. Neglecting Cabibbo-suppressed and penguin-induced transitions, such 
mechanisms are: 

• the 6-quark transition b c , with the c quark having the role of spectator; the 
corresponding final states are of the kind J/ip n, J/ip iu; 

• the charm quark transition c — > s W~, with b as spectator and possible final states 
Bg n, Bg tv, etc.; 

• the annihilation modes cb — W~ . 

The first two mechanisms are responsible of the largest part of the B^. decay width 
P, 1^, 0. In particular, the measurement (^ provides us with an indication that the 
dominant B^. decay mechanism is the c-quark decay, which implies a B^. lifetime in the 
range tb, = (0.4 — 0.7) ps 0, whereas dominance of the 6— quark decay mechanism would 
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produce a longer lifetime: tb^ = (1.1 — 1.2) [^]. Various analyses of Be transitions 
induced by the two mechanisms are available in the literature for example, a QCD 
sum rule calculation of semileptonic decays suggested the dominance of the charm 
transition p. 

As far as the annihilation processes are concerned, the leptonic radiative decay mode 
Be —>■ and the leptonic decay without photon in the final state represent a minor 
fraction of the Be full width. Nevertheless, their analysis is of particular interest, both 
from the phenomenological and the theoretical point of view. From the phenomenological 
side. Be annihilation modes are governed by Vet] therefore they are Cabibbo-enhanced 
with respect to the analogous Bu decays, and represent new channels to access this matrix 
element. From the theoretical viewpoint, the purely leptonic and the radiative leptonic 
Be transitions are interesting since, in the nonrelativistic limit of the quark dynamics, 
both their rates can be expressed in terms of a single hadronic parameter, the Be leptonic 
decay constant /^^ P, |T0[. In this limit, a relation between the widths of the processes 
Be — > £z/7 and Be (.v can be worked out p, p!0|] : 



a m 
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with ri = Eq.(0) implies that the width of the radiative leptonic Be decay into 

4:71 mf 

muons is nearly equal to the purely leptonic width T{Be fiiy), whereas, in the case of 
electrons in the final state, the radiative leptonic mode is largely dominant. 

Eq.(^) presents uncertainties coming from the used values of the charm and beauty 
quark masses. Moreover, there could be corrections if the quark dynamics in the Be 
meson deviates from the nonrelativistic regime, and the size of the corrections is useful 
for understanding the theoretical uncertainty affecting the ratio (|^). 

Corrections to the ratio can be estimated by considering a model for the Be meson 
where relativistic effects in the constituent quark dynamics are, at least partially, taken 
into account; the present letter is devoted to such a study. 

In order to analyze the decay mode Be —>■ {i = e, fi), we follow the method adopted 



in ref.[|TT[ to investigate the analogous B^ transition. Let us consider the process 

B;ip)^rip,)uip2)7{k,e) (4) 
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whose amplitude can be written as 

A{B, ^ £z/7) = e*' {L>^ ■ U,, - lefByL,,) , (5) 

with Gp the Fermi constant and e the photon polarization vector. The current L^^ = 
i{pi)'~f^{l — 75)2^(^2) represents the weak leptonic current in (^). The hadronic function 
n^^(p, k) is the correlator 

Tif.Ap, k)=t J < 0|T[J^(x)K(0)]|5,(p) > (6) 

with J^{x) = c(x)7^(l — 75)6(x) the weak hadronic current governing and given by 
Vu{0) = e[Qc c(0)7i/c(0) +Qb &(0)7i,6(0)], eQc and eQb being the charm and beauty quark 
electric charges. The momentum q is defined as q = pi + p2- Therefore, the first term in 
d^) correspond to the photon emitted from the meson, and the second term represents the 
contribution of the photon emitted by the charged lepton leg. Notice that the leptonic 
constant Jb^ is defined by the matrix element 

< 1 07^756 1 5, (p) >= ifB^p^, (7) 



and that L^^, reads 



1 

:^^ + ^ - mi 



V = ^ipih" . , . ^. 7^(1 - 75)^(^2) 



The hadronic function Ilfj,u{p, k) can be expanded in independent Lorentz structures 

^iiu{p ■k) = a Pf,p^ + (3 kf,k^ + ( k^p^ + 5 p^k^ + ^ Qt^u + ive^^p^pfk"" (9) 

(with the invariant functions a, . . .rj depending on p- k) so that the requirement of gauge 
invariance for the amplitude implies the condition: 

{a + O{p-k)+C-tefB^ = . (10) 

We shall see in the following that this condition is satisfied in our calculation. 
The final expression of the amplitude, eq.(^, 

A{B, ^ M) = e*'' L^[{a + () {k^p, -p-k g^,) + i r] e^^^^p^k'^] e*^ , (11) 

is given in terms of the form factors r] and C + c^^ related to the vector and the axial vector 
weak current contributions to the process (Hf), respectively. 



In order to compute the invariant functions in (|TTp we employ the costituent quark 
model developed in ref.|T^ to describe the static properties of mesons containing heavy 
quarks. In the correlator we write the state of the pseudoscalar [be] meson at rest in 
terms of a wave-function ips^ and of quark and antiquark creation operators: 

\B;>=^^-^^Jdk^^BMbKk,r,a)c^-k,s,P)\0> ; (12) 

a and (3 are colour indices, r and s spin indices; the operator creates a b quark with 
momentum k, while creates a charm antiquark with momentum —k. The wave-function 
ipBcy describing the quark momentum distribution in the meson, is obtained as a solution 
of the Salpeter equation 



yk^ + ml + ^JP + ml-mB^}4jBAk) + J dk' V{k,k') ^jB^k') = (13) 

stemming from the quark-antiquark Bethe-Salpeter equation, in the approximation of an 
istantaneous interaction represented by the interquark potential V. Within the model in 
is chosen as the Richardson potential [O, which reads in the r— space: 



V{r)= A[Ar-^] , (14) 



33 - 2n 



f 



Ar 



with A a parameter, nf the number of active flavours, and the function f{t) given by: 

^ 4 /■°° , sin(qt) r 1 1 1 ^, 

71 Jo q ^ln{l + q'^) g^J 

The linear increase at large distances of the potential in (|T^ provides QCD confinement; at 
short distances the potential behaves as — with as{r) logarithmically decreasing with 
the distance r, thus reproducing the asymptotic freedom property of QCD. A smearing 
procedure at short-distances is also adopted, to account for effects of the relativistic 



kinematics [0. Finally, spin interaction effects are neglected since in the case of heavy 



mesons the chromomagnetic coupling is of order of the inverse heavy quark masses. 

The interest for eq.(|T3|) is that relativistic effects are taken into account at least in 
the quark kinematics. Therefore, the same wave-equation can be used to study heavy- 
light as well as heavy-heavy quark systems, and one case can be obtained from the other 
one by continuously varying a parameter, the u{c) quark mass. From the analysis of the 
solutions, a comparison between the two cases, the heavy-light and the heavy-heavy quark 
meson systems, can be meaningfully performed. 



4 



Eq.(|TB|) can be solved by numerical methods, as described in |T^, fixing the masses 
rric and rrib of the constituent quarks, together with the parameter A, in such a way that 
the charmonium and bottomomium spectra are reproduced. The chosen values for the 
parameters are: rrib = 4.89 GeV and = 1.452 GeV, with A = 397 MeV pI. A fit of 
the heavy-light meson masses also fixes the values of the constituent light-quark masses: 
TTLu = rrid = 38 MeV [Q. For the be system all the input parameters required in (|13D 
are fixed from the analysis of other channels, and the predictions do not depend on new 
external quantities. The numerical solution of (p!3D produces spectrum of the be bound 
states; the predicted masses of the first three radial S*— wave resonances are reported in 
Table 1. 



radial number 


mass (GeV) 


leptonic constant (MeV) 





6.28 


432 


1 


6.80 


313 


2 


7.16 


265 



Table 1: Spectrum and leptonic constants of the pseudoscalar be mesons. 



The spectrum in Table 1 agrees with other theoretical determinations based on con- 



stituent quark models |IJ]. As for QCD sum rule and lattice QCD results, in the value 
ttib^ = 6.35 GeV was found using two-point function QCD sum rules, whereas a recent 
lattice QCD analysis predicts: m^^ = 6.388 ± 9 ± 98 ± 15 GeV, with the larger error 



related to the quenched approximation ||T^ . Within the errors, the mass of the Be meson 
in Table 1 agrees with the CDF measurement reported in eq.([^). 

Also the wave-function ipBc can be obtained by solving eq . (p!3D . We use the covariant 
normalization: 

^ J dkl^PsM' = '^rriB^ (16) 
and define, in the meson rest-frame, the reduced wave-function UB^{k) {k = \k\): 

usXk) = (17) 

which is normalized as: \^ dk\uBS}^^^ = '2mB^- In fig.l the Be meson wave-function 
UB^k) is depicted, together with the function ub^ of the Bu meson. 
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Figure 1: Be meson wave- function (continuous line). The wave-function of the Bu meson 
is also depicted (dashed line). 



The meson wave-functions, together with the predictions of the spectrum of the bound 
states, are the main results of the model; they allow us to calculate hadronic quantities 



such as the leptonic decay constants and the strong couplings to light mesons |12, 16 
In particular, using UB^{k) depicted in fig.l, we can infer the size of the deviation from 
the nonrelativistic limit in the Salpeter equation (|TB|). As a matter of fact, the average 
squared quark momentum < k"^ > turns out to be < A;^ >= 0.95 GeV^, and the ratios 
< k"^ > /m^ = 0.43 and < k"^ > /ml = 0.04. In the case of the Bu meson, the average 
squared quark momentum < fc^ > is < /c^ >= 0.47 GeV^, to be compared to = 
1.4 X 10^^ GeV^. Therefore, in the case of Be mesons, deviations from the nonrelativistic 
limit, although small, are not negligible, mainly due to term related to the charm quark. 

Before coming to the analysis of the radiative leptonic Be decay and to the calculation 
of the ratio (|^), let us evaluate the leptonic constants /^n of the 5- wave excitations, 
defined by matrix elements analogous to (|^. The numerical results, obtained from the 
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eigenfunctions of eq.(|IB|), are reported in Table 1. Few comments are in order. The value 
of in Table 1 is compatible, within the theoretical errors, with the result Jb^ = 360 ±60 
MeV obtained from QCD sum rules (the uncertainty related to the input parameters 
of the potential model is estimated of ±15% for the leptonic constants |ll2[)- On the 



other hand, the leptonic constants in Table 1 are smaller than the outcome of several 



constituent quark models considered in [|14|, based on a purely nonrelativistic description 
of the be system. Notice that the decreasing values of /^n, for increasing radial number 
n, is due to the nodes of the wave-functions of the radial excitations. Finally, the Be 



leptonic constant turns out to be compatible with a lattice NRQCD determination [jT7|. 



Let us now consider the correlator (P). Writing the state |-Bc > as in (|T^) and expand- 
ing the T-product according to the Wick theorem, we can express the product of the two 
currents J^(a;)V^(0) in terms of quark creation and annihilation operators; then, exploiting 
the anticommutation relations between such operators, we can derive an expression for 
n^jy in terms of the Be wave function, ub^, which is analogous to the expression reported 
in [|ll|] for the B^ decay. 

It is useful to relate the invariant functions in (j^) to the various components of the 
hadronic tenson II^j,. In the B^ rest-frame p = {mB^,0) and choosing k'^ = (A;°, 0, 0, A;°) 
one gets from eq.(^: 



V{k') = -^^^,Ill2{k') . (18) 
Therefore, the condition (|I0|) ensuring gauge invariance can be written as 

^e/B. = -^(noo±n33-n3o-no3)±n3o-n33 , (i9) 

rriB, 

a condition that must be checked in our analysis. The explicit calculation, using the 
expression of II^i, in terms of the meson wave-function shows that eq. (|l9|) is verified 
provided that the leptonic constant /^^ is given by 

/bc = V3- / dk' k' UB,.{k')[^ -P- ^ p l-TiT— wTT— r 20 
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with Eb^c = yk''^ + itlI^. Indeed, eq.(^) is the expression for fs^ obtained in the frame- 
work of the constituent quark model |T^, hence the gauge invariance property of the 
amphtude (|) is preserved in our calculation. f\ 

Having checked the property of gauge invariance, we can simply compute the decay 
width of the mode (H) and the photon energy distribution. Notice that we compute the 
photon energy spectrum for a photon energy larger than 100 MeV, which represents the 
photon energy resolution we assume in our analysis. 

The expression for the width of the decay (H) , considering massless leptons in the final 
state {me ^ 0), is 

(72 IT/, |2 /•mBr/2 

r(i?e ^ M) = dk\\mB^ - 2k')[\Un + tefnf + lUul'] . (21) 

Using the data in (|]),(|D, together with Vcb = 0.04, we predict: 

r(fic M) = 4 IQ-^^GeV , B{B, iu-f) = 3 10^^ . (22) 
The energy spectrum of the photon, computed using eq.(pl|), is depicted in fig. 2; it looks 



symmetric with respect to the point = as observed also in |T0|, |T9 



The result (^) confirms that the rate of the radiative leptonic Be decays into muons 
(and electrons) is sizeable, and could be accessed at the future hadronic machines, such 
as LHC. 

Let us now consider the ratio (§). In order to determine it, we need the width of the 
purely leptonic Be decay into electrons and muons, given by 



r(B, ^ fo) = ^\V^pfi(j^Y,n%J 1 - . (23) 



which depends on fs^- Using the value in Table 1 we obtain: 



T{Be ^ /iz/) = 1.1 10-^^ GeV , B^B, /iz/) = 7.8 10"^ , (24) 

corresponding to Rn = — t~ ^ — This result must be compared to the value 

T{Be /iz/) 

~ 0.8 obtained using eq.(|^) ||, 0], and suggests that the corrections to the nonrela- 
tivistic limit, although small, play a role in modifying the prediction (H). In any case, the 



■^In the case of the radiative leptonic i?„ decay, the contribution proportional to /b turns out to be 
numerically negligible. 
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Figure 2: Photon energy spectrum in the decay 



widths of the radiative leptonic and the purely leptonic Be transitions are still compara- 
ble, a result which is very different from the case of radiative leptonic Bu decays, which 
are enhanced with respect to the purely leptonic decay modes into muons and electrons. 
In that latter case, the helicity suppression, displayed by the factor {^^)^ in (^) and 
producing very small decay rates, is efficiently avoided by the presence of the photon in 
the final state also in the case of muons, so that the branching fraction of the radiative 
leptonic B^ modes into muons turns out to be enhanced by one order of magnitude with 
respect to the purely muon leptonic decay |TT|, |2D|. 

As for the decay into electrons, we predict B{Be —>■ eu^) ^ 2 x 10^^. 

The decay mode Be has been investigated, as mentioned above, in the frame- 

work of the nonrelativistic quark model (NRQM) |P, |lOl, and also using the light front 
quark model (LFQM) p| and light-cone QCD sum rules (LCSR) pf. In Table 2 we 
report the various predictions, obtained using the value of fs^ reported in Table 1. We 
also report the estimates of the ratio 



(25) 
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which represents the relative fraction of the final state £1/7 coming from the different 
sources Bu and Be- Our small result: r = 0.03 means that a large number of iwy final 





r 


method 


ref. 


3 X 10"^ 


0.03 


Rel. QM 


This paper 


6 X IQ-^ 




NRQCD 


|10| 


6 X 10-^ 


0.06-0.1 


NRQM 


i 


2 X 10-5 




LFQM 


|18| 


2 X 10^5 


0.1 


LCSR 


|19| 



Table 2: Results for B{Bc — > iw-f) by different approaches, 
states should come from Be decays; the actual fraction of radiative leptonic final states can 



be obtained by multiplying eq. (|25D by the probabilities of producing Bu and Be mesons 
from b quarks. 

Let us conclude our analysis, based on a QCD-inspired relativistic constituent quark 
model, observing that the quark dynamics inside the Be meson could modify the pre- 
diction (^). Therefore, the range [0.4 — 0.8] for the ratio (§) can be interpreted as the 
theoretical uncertainty for this quantity. The measurement of the radiative leptonic Be 
decay rate, together with the purely leptonic rate, although challenging from the experi- 
mental viewpoint, is one of the expected results in the LHC era. 



Acknowledgment s 

We thank G. Nardulli and N. Paver for discussions. 



10 



References 



[1] 

[2] 



CDF CoUab., F. Abe et al, Phys. Rev. Lett. 81 (1998) 2432; Phys. Rev. D 58 (199^ 
112004. 



For a review of the properties and decays of Be see: S.S. Gershtein et al., [hep- 
ph/9803433| , |hep-ph/95043T9 . 



[3] K. Cheung, Phys. Rev. Lett. 71 (1993) 3413. 

[4] K. Kolodziej, A. Leike and R. Riickl, Phys. Lett. B 355 (1995) 337 C.H. Chang, 
Y.Q. Chen, CP. Han and H.T. Jiang, Phys. Lett. B 364 (1995) 78; C.H. Chang, 
Y.Q. Chen, and R.J. Oakes, Phys. Rev. D 54 (1996) 4344.; A.V. Berezhnoy, V.V. 
Kiselev, A.K. Likhoded and A.I. Onishchenko, Yad. Fiz. 60 (1997) 1886. 

[5] M. Lusignoh and M. Masetti, Z. Phys. C 51 (1991) 549. 

[6] C. Quigg, in Proceedings of the Workshop on B Physics at Hadron Accelerators, 
Snowmass, Colorado 1993, edited by P. McBride and C.S. Mishra, pag.439. 

[7] M. Beneke and G. Buchalla, Phys. Rev. D 53 (1996) 4991. 

[8] P. Colangelo, G. Nardulh and N. Paver, Z. Phys. C 57 (1993) 43. 

[9] C.H. Chang, J.P. Cheng and CD. Lii, Phys. Lett. B 425 (1998) 166. 

[10] G. Chiladze, A.F. Falk and A.A. Petrov, Phys. Rev. D 60 (1999) 034011. 

[11] P. Colangelo, F. De Fazio and G. Nardulh, Phys. Lett. B 386 (1996) 328. 

[12] P. Cea, P. Colangelo, L. Cosmai and G. Nardulh, Phys. Lett. B 206 (1988) 691; P. 
Colangelo, G. Nardulli and M. Pietroni, Phys. Rev. D 43 (1991) 3002. 

[13] J. L. Richardson, Phys. Lett. B 82 (1979) 272. 

[14] W. Kwong and J. Rosner, Phys. Rev. D 44 (1991) 212; E. Eichten and C. Quigg, 
Phys. Rev. D 49 (1994) 5845; S.S. Gershtein et al., Phys. Rev. D 51 (1995) 3613; 
L.P. Fulcher, Phys. Rev. D 60 (1999) 074006 and references therein. 

[15] H.P. Shanahan, P.Boyle, C.T.H. Davies and H. Newton, UKQCD Collab., Phys. Lett. 
B 453 (1999) 289. 

11 



[16] P. Colangelo, F. De Fazio and G. NarduUi, Phys. Lett. B 334 (1994) 175. 

[17] B.D. Jones and R.M. Woloshyn, Phys. Rev. D 60 (1999) 014502. 

[18] C.C. Lih, C.Q. Ceng and W.M. Zhang, Phys. Rev. D 59, 114002 (1999). 

[19] T.M. Ahev and M. Savci, Phys. Lett. B 434 (1998) 358; J. Phys. G 25 (1999) 1205. 

[20] G. Burdman, T. Goldman and D. Wyler, Phys. Rev. D 51 (1995) 111; D. Atwood, 
G. Eilam and A. Soni, Mod. Phys. Lett. A 11 (1996) 1061; G.Eilam, I. Halperin 
and R.R.Mendel, Phys. Lett. B 361 (1995) 137; P. Colangelo, F. De Fazio and G. 
NarduUi, Phys. Lett. B 372 (1996) 331. 



12 



